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ABSTRACT 
A high resolution ocean general circulation model with a horizontal resolution of 10 Km is used to investigate salinity in the 
Northern Indian Ocean [65E-95E; 5N-20N]. The model simulated salinity during the years 1998-2014 is compared against the 
satellite and in-situ observations. The model is able to capture the contrasting distributions of sea surface salinity (SSS) in the 
Arabian Sea and the Bay of Bengal. The Bay of Bengal salinity values are generally lower than the Arabian Sea. The model is also able 
to reproduce the strong seasonality of the sea surface salinity in the changing climate. The inter-annual SSS variability is 2-4 times 
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smaller than seasonal variability in the region of study. The regions of large inter-annual variability are located near river mouths in 


the Bay of Bengal. The inter-annual variability of the SSS in the northern Indian ocean is also examined in the light of recent changes 
in the El Nino-Southern Oscillation (ENSO), Indian Ocean Dipole (IOD) during the period of study. The relative role of wind and solar 
radiation on the sea surface salinity is also investigated. 
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1. INTRODUCTION 

The northern Indian Ocean is forced by semi-annually reversing winds that blow from the south west during the summer monsoon 
and from the north-east during the winter monsoon. The Strong evaporation over precipitation and inflow from the Red Sea and 
Persian Gulf makes the Arabian Sea (AS) much saltier (Antonov et al., 2010; Chatterjee et al., 2012). In contrast to this, the Bay of 
Bengal (BoB) experiences high precipitation and river discharge by adjoining river Ganga, Brahmaputra and Irrawaddy during the 
monsoon and post-monsoon months (Rao et al., 2003; Sengupta et al., 2006; Durand et al., 2011). The river runoff is entering initially 
through mouth of river and advected into the southern BoB and then in the eastern AS by the East India Coastal Current 
(Vinayachandran et al., 2007). The yearly freshwater flux (precipitation plus river runoff) received by the BoB is largely exceeds the 
freshwater flux evaporated back to the atmosphere (Shenoi et al., 2002; Sengupta et al., 2006). As a result, sea surface salinity (SSS) 
of northern BoB is very low (Parampil et al., 2010; Han et al., 2001) and highly stratified. These fresh near surface waters affects the 
depth of wind-induced vertical mixing which is responsible for shallow mixed layers in the northern BoB (Lukas and Lindstorm 1991, 
de Boyer Montegut et al., 2007, Girishkumar et al., 2013). Due to high stratification in northern BoB, formation of barrier layer takes 
place which prevents the entrainment of cool, deep water in to surface which is responsible for making warmer upper ocean 
temperature of BoB. The salinity stratification helps in reducing the cooling effects during cyclone in the BoB, which prone to 
intensification of tropical cyclones (Sengupta et al., 2008; Neetu et al., 2012). This low SSS and stratification in the BoB is also 
influence the regional climate. It also affects the climate of the Indian subcontinent including the summer monsoon rainfall thus 
affecting lives of billions of people living there. Salinity also modulates the intra-seasonal variability of the SST in the Northern Bay 
of Bengal (Vinayachandran et al., 2012).The SSS observation in Northern Indian Ocean is very sparse in space as well as in time. For 
better description of the SSS seasonal evolution within the Northern Indian Ocean, and to assess all the processes that govern the 
seasonal SSS evolution, the observation data are not yet quantitatively sufficient. So Numerical ocean circulation models have the 
distinct advantage of generating long time series of salinity at high spatial and temporal resolution which can help in understanding 
processes responsible for salinity variability. The role of river runoff, precipitation, and advection of fresh water from northern BoB to 
AS through EICC have been identified through sensitivity experiment using modelling study (Diansky et al., 2006; Wu et al., 2007; 
Sharma et al., 2010; Girishkumar et al., 2013; Benshila et al., 2013; Akhil et al., 2014; Chaitanya et al., 2015). The seasonality of sea 
surface variability is strongly influence by monsoonal current (Shankar et al., 2002) which advect low salinity water of BoB into AS 
during Winter monsoon season and high salinity water of AS into BoB during summer monsoon season. The salinity is important 
parameter for understanding the ocean hydrological cycle, which is least understood and important component of the climate 
system (Schmitt, 1994), it plays an important role in controlling near-surface vertical mixing and affect several climatic phenomena, 
i.e. Tropical cyclone, monsoon rainfall (Sanilkumar et al., 1994; Vinayachandran et al., 2015). 

In this paper, we describe the surface and subsurface salinity variability on seasonal and interannual time scales using a high 
resolution ocean model and estimate the contribution of wind and solar radiation in modulation of the salinity. In section 2, we 
describe the configuration of the ocean model used for the study and details on setup of sensitivity experiments. In section 3, we 
describe the validation of model results. In section 4, we describe the surface and subsurface salinity variability on the seasonal and 
interannual time scales. We describe the contribution of wind and solar radiation on salinity using the sensitivity experiments in 
Section 5. Conclusions are in Section 6. 


2. MODEL CONFIGURATION 

We use the MIT General Circulation Model (MITgcm) (Marshall et al., 1997) to carry out the present study. The MITgcm is a z- 
coordinate model, which uses hydrostatic approximation. The model space is discretized using staggered Arakawa C grid. We 
configure the model at a high horizontal resolution of 10 km in the Bay of Bengal region around [65E-95E; 5N-20N] during 1 Jan 
1998 to 31 Dec 2014. The highest vertical resolution is taken as 5 m near the ocean surface and then gradually telescoping out at 
depth. We used a total of 28 vertical levels with maximum vertical spacing of 500 m near the ocean floor. The model uses open 
boundary conditions on all four sides, which are taken from the ORAS4. We used the KPP non-local vertical mixing scheme (Large et 
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al., 1994) to resolve the sub-grid scale processes. The eddy (harmonic) viscosity and diffusivity satisfying the CFL stability criterion is 
chosen in horizontal and vertical. We use no-slip condition on sides as well as bottom. The bottom frictional drag coefficient is taken 


as 0.001. A 3rd order direct space-time advection scheme is employed for temperature and salinity. The bathymetry of the region of 
interest in the BoB (Figure 1) is taken from the Smith and Sandell (1997) 1’ bathymetry data. The nonlinear equation of state is used 
following the Jackett and McDougall (1995). The initial temperature and salinity are derived from the World Ocean Atlas 2013 
(WOA13) data (Locarnini et al., 2013; Zweng et al., 2013). The 6-hourly NCEP/NCAR reanalysis (Kalnay et al., 1996) forcing, namely, 
2m air temperature, humidity, downward long-wave and short-wave radiations, precipitation, zonal and meridional winds and runoff 
are used to force the model. Evaporation are calculated by model itself using bulk formula. Model is relaxed with monthly 
climatology of SST and SSS from WOA13. 


Figure 1 Bathymetry of the region of study 


3. QUALITY OF MODEL OUTPUT 

We begin our analysis with evaluation of model SSS against the North Indian Ocean Atlas (NIOA) climatology provided by Chatterjee 
et al. [2012]. Figure 2 shows the climatology and standard deviation of SSS from model during 1998-2014. The corresponding maps 
from the NIOA Climatology are also shown. It is clear from the figure that model is able to capture the contrasting structure of the 
sea surface salinity (SSS) in the AS and BoB. Low SSS in the BoB and high SSS in the AS is primarily due to excessive evaporation over 
precipitation. In the BoB, the lowest SSS is observed in the northern domain, which gradually increases towards south. In the AS, 
high SSS is found in northern domain which decreasing values towards south. The strong SSS variability (>1 psu) is found in 
northern BoB especially near the river mouth and along the east coast of India. The south east Arabian Sea region is another 
relatively high SSS variability region. 
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Figure 2 Mean (1998-2014) and standard deviation of the monthly climatology of SSS from model and NIOA 
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Figure 3 shows comparison of seasonal climatology obtained from model with NIOA Climatology (Chatterjee et al. 2012). It is 


clear from the figure that model performs reasonably well in simulating the seasonal SSS variability of the region of interest. The 
northern BoB is a region marked with high freshwater flux during the summer monsoon season (June-July-Aug; JJA). The freshening 
extends southward and becomes maximum during the post monsoon (Sept-Oct-Nov; SON) months and as a result the SON seasons 
is also showing the lowest SSS in the BoB due to of high runoff from major Indian rivers. 
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Figure 3(a) Seasonal climatology obtained from model and NIOA for the year 1998-2014 
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Figure 3(b) Seasonal climatology obtained from model and Aquarius Satellite for the year 2011-2014 


During winter (Dec-Jan-Feb; DJF) season, the freshwater moves further south along the east coast of India through the EICC and 
gets advected in the AS through westward-flowing Winter Monsoon Current (WMC) (Jensen, 2001, Shankar et al., 2002, McCreary 
2001; Sengupta et al. 2006). This advection of warm low-salinity water over cool high-salinity water (Shankar et al., 2004; Durand et 
al., 2004; Gopalakrishna et al., 2005) is responsible for formation of temperature inversions in south east AS (Kurian et al.2006). The 
salty water from the AS is advected into the BoB during the summer monsoon months (Murty et al.1992; Vinayachandran et al. 
1999). The highest SSS in the AS is found during the MAM. This is due to much higher evaporation than precipitation during this 
season. These seasonal features of SSS are well captured by the model and also in conformity with the past studies (Rao et al., 2003, 
Chaitanya et al. 2014). The high-resolution model also captures the small-scale structures such as the narrow winter monsoon 
freshwater strip along the EICC. However, this known feature is absent in the NIOA. This may be due to coarse resolution of the 
NIOA and associated interpolation errors (Chatterjee et al. 2012, Chaitanya et al. 2014). 

We also compare model SSS with Aquarius satellite data for the available period. Figure 3(b) shows seasonal climatology 
obtained from model as well as from Aquarius Satellite for the year 2011-2014. From the figure it is clear that model is able to 
reproduce the satellite SSS variability. In order to see the model's performance in getting the correct subsurface salinity, we compare 
the model output with NIOA. Figure 4 shows the depth-latitude plot of mean salinity during 1998-2014 at 70E and 90E longitudes. It 
is clear from figure that model is able to capture the layered structure of sub-surface salinity. 


4. SEASONAL VARIABILITY OF SSS 
Figure 5 shows the seasonal SSS variability (standard deviation) over the AS and BoB regions. The high seasonal variability is seen in 
the northern BoB and eastern coast of India, especially near the major river mouths. In AS, highest variability is found in the southern 
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AS during DJF and MaM. During SON, the entire AS shows very low SSS variability. However, it is clear from the figure that northern 
AS shows relatively small salinity variability as compared to the southern AS in all seasons. In MAM, the southern BoB also shows 
high SSS variability. Similarly high SSS standard deviation is also seen along the southeastern coast of India during SON. The high 
SSS variability observed in the northern BoB in the summer monsoon months of JJA, and post-monsoon SON months are due to 
increased freshwater flux (river runoff) from major adjoining rivers of the region, namely, Ganga, Brahmaputra and Irrawaddy. 
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Figure 6 partial correlation of SSS with different climate index 


paoe | vas 


ARTICLE 


5. RELATION OF SSS WITH DIFFERENT CLIMATE MODES 
In order to examine variability of the SSS patterns over the AS and BoB vis-a-vis variability of the known climate indices, we find the 
spatial correlation (Figure 6) of SSS with the time series of NINO3.4, Dipole Mode Index (DMI), Southern Dipole Index (SDI) and 
Indian Monsoon Index (IMI). 

From figure 6, it is clear that SSS variability over the AS and BoB is influenced by the ENSO, IOD, SIOD and monsoon variability. 
The correlation of NINO3.4 is highest in the AS in all seasons. DMI has highest correlation in autumn in eastern coast of India and 


eastern Arabian Sea with correlation ranging from 0.6 to 0.75. On the other hand, the south eastern BoB and south western AS show 
negative DMI correlation in the range of -0.6 to -0.9. 


6. CONCLUSION 

The goal of present study is threefold: (i) to demonstrate the quality of salinity derived from the high resolution general circulation 
model in the Arabian Sea and Bay of Bengal; (ii) to explore the variability of surface and subsurface salinity; (iii) to assess the role of 
near surface wind and solar radiation on the seasonal variability of SSS. The SSS from model is found to be in good agreement with 
the satellite and in-situ observations. The model is able to reproduce the distinct SSS features in the AS and BoB. The model also 
captures small-scale structures such as the narrow winter monsoon freshwater strip along the EICC. This feature is seen only in the 
high-resolution model output, whereas, it is absent in the coarse resolution reanalysis data. The inter-annual SSS variability is 2-4 
times smaller than seasonal variability in entire study domain. The regions of large SSS variability are located near the river mouths 
in the Bay of Bengal. The SSS also gets influenced by different climate modes. The spatial correlation of SSS with NINO3.4 is higher 
in AS than BoB for all seasons. The DMI has highest correlation in autumn along the eastern coast of India and eastern Arabian Sea. 
On the other hand, south eastern BoB and south western AS show negative correlation with SSS. The SSS is also strongly affected by 
the wind and solar radiation. The SSS decreases in the absence of wind as well as solar radiation. 
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